Resonant magneto-conductance through a vibrating nanotube 
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We address the electronic resonant transport in presence of a transverse magnetic field through the 
single level of a suspended carbon nanotube acting as a quantum oscillator. We predict a negative 
magneto-conductance with a magnetic-field induced narrowing of the resonance line and a reduction 
of the conductance peak when the nanotube is asymmetrically contacted to the leads. At finite bias 
voltage we study the threshold for phonon-assisted transport. 

PACS numbers: 73.63.-b, 71.38.-k, 73.63.Kv, 85.85. 



A distinctive non-classical feature of any quantum 
state is the possibility that part of the system is spa- 
tially delocalized. Probing delocalization in macroscopic 
systems is important to validate quantum mechanics on 
that scale. For large molecules, this has been done by the 
observation of interference fringes with diffraction exper- 
iments Is it possible to observe quantum delocaliza- 
tion for a mechanical oscillator? Recently, it has been 
predicted that the magneto-conductance can be used as 
a detector of the quantum delocalization of a suspended 
carbon nanotube due to the Aharonov-Bohm effect on 
the electrons crossing the device [Ij. This spectacular 
prediction was done in the tunnel regime. Though sim- 
ple and transparent from the technical point of view, this 
regime is not optimal for the experimental observation of 
the effect for two reasons: (i) the current is very low, (ii) 
electrons can interfere only once since a single crossing 
through the device contributes to the current. 

In this article, we consider the resonant transport 
through a single electron state of the nanotube. At 
resonance, the electron channel is fully open for a 
static nanotube. One could expect that the magneto- 
conductance signal for the suspended nanotube will be 
greatly enhanced: The electrons bounce many times in- 
side the structure before leaving, therefore allowing mul- 
tiple interference. Thus, even if the phase acquired at 
each passage is small, the accumulated phase can be 
large. By performing a calculation with Keldysh non- 
equilibrium Green's function technique at lowest order 
in the electron-phonon coupling, we find the following 
results: The shape of the resonance as a function of the 
gate voltage is modified by the magnetic field. At reso- 
nance and for vanishing temperature, the linear conduc- 
tance depends on the magnetic field only if the coupling 
to the leads is asymmetric, while the resonance width 
is reduced by the magnetic field. The current-voltage 
characteristics shows a magnetic-field-dependent singu- 
larity at the threshold of one-phonon absorption. These 
prominent features constitute a measurable signature of 
the quantum delocalization of the nanotube. 



We model the system with the following Hamiltonian: 

= E E ^kA.^k. + £d + nwhh + HT. (1) 



u—l^r k 

Here, a^^ (<^j.r) ^ creation operator for the electronic 
single particle states in the left (right) lead. The en- 
ergy spectrum in each lead is ^kv = £k — i^u, where the 
difference of the chemical potentials /i^ — /i; = eV is re- 
lated to the bias voltage V. The leads are connected 
by a suspended nanotube placed in a strong magnetic 
field perpendicular to the nanotube's oscillation plane. 
We single out the fundamental bending mode with reso- 
nance frequency lo for which is the creation operator of 
quantum excitations. The single relevant electronic level 
in the nanotube for which is a creation operator sits 
at energy £d that can be controlled with an external gate 
voltage. For simplicity, we consider fully spin-polarized 
electrons (having in mind the large magnetic field for 
which our theory applies) and neglect a possible orbital 
degeneracy in the nanotube. The last term in Eq. ([T|) 
models electron transfer from the leads to the nanotube 
in presence of the magnetic field 01 : 



Hj 



h.c. 



(2) 



Here, a; = 5 + 6^ is the displacement operator of the nan- 
otube in units of the amplitude uq for zero-point fluctua- 
tions, ti and tr are the tunneling matrix elements at the 
point contacts between the nanotube and the leads. The 
factor a = ai = — = gBLuo/^o is the magnetic flux 
(in units of the flux quantum <f>o = h/e) through the area 
swept by the nanotube with length L in the ground-state 
quantum fluctuation at magnetic flcld B {g is a numer- 
ical factor of order one). Typically, a is very small: for 
a single wall nanotube of length L = l/xm, (ci;/27r) = 500 
MHz, one flnds mq ~ 1 pm. With B = 40T, one ob- 
tains a ~ 0.1. The Hamiltonian ([T]) resembles the one 
proposed to study polaronic transport through a vibrat- 
ing molecule [3, [j| (there ai = a^, £d includes a pola- 



2 



ronic shift — a^w, and x stands for dimensionless mo- 
mentum), but they lead to qualitative different behavior. 
The Hamiltonian of Ref. [5| describing Coulomb block- 
ade at resonant tunneling reduces to ([T]) in some limit. 
We address a regime that was not considered there. 
The current operator in the left lead is / = 



ite/h)[H,Ni], with Ni = Efc 



In the stationary 



regime, the dc current flowing through the device is 

^=-f E(*^(e^"'«U)-c.c.) , (3) 

k 

where the brackets denote a quantum-statistical aver- 
age. In order to evaluate the current, we use the 
Keldysh theory for non-equilibrium systems 0]. We 
define the retarded, advanced, and Keldysh electronic 
Green's functions: G^/„^(t) = TiO{±t){{an{t),aliQ)}), 
Gn,n'it) = -i{[an{t),aliO)]) (with n,n' = kl,k'r,d) 
and we build a triangular 2X2 matrix G out of them, 
with elements Gn = G^, G12 = G'^ , G22 = G^, 
and G21 = 0. We define similarly a Green's func- 
tion Gxd.ki related to Eq. such that, for instance, 

GSt(0 = T*0(±t)({arf(t)e--(*),aL(O)}). The rela- 
tion Gki,k'i{£) = gki{£)Sk,k' + gki{£)tiGxd,k'i{£) holds in 
Fourier space, where gki is the Green's function in the un- 
coupled left lead (at H-r ~ 0). We introduce = J^k d^v 
and Gu^^' — "^f^ t^i Gkv,k'v' ■ In the wide-band limit: 

gu^^ = TiT'Pi^, g^y = 2[1 - 2ni]g^ , where are the 
densities of states in the leads and = np{e — /i^), 

with np the Fermi distribution function at temperature 
T. Then, Eq. ^ can be rewritten: 



e 

-Re 
h 



gT\e)Gi,i{e) 



(4) 



In the absence of electron-phonon coupling (at a = 0), 
the Green's function on the dot is known: G^''|/*^'^(e) = 



{e - ed±iT)-\ while G|;°|f (e) 
r']~'E.=/,.r,[l - 2n.(£)]. Here, 
Tf/r = T^Pi/r\ti/r\'^ givc the broadening of the resonant 



= ~2i[{e ~ Edf + 
F = + F,. and 



level due to its hybridization with the leads. Then, one 
gitiG^^\t*igi. Inserting this into Eq. (g]). 



gets: Gf} 



gi 



de[ni{e) - nr{e)]T{e), 



(5) 



one retrieves the result 
h 

with the elastic Breit-Wigner transmission coefficient 
through the non-interacting resonant level 

4F,F, 



T{e) = 



F2 



(6) 



In particular, the linear conductance at resonance and 
T = 0, Gmax = (e^/ft.)4F;Fr/F2, reaches the conductance 
quantum for symmetric contacts (F; = F^). 

We now consider the coupling with phonons perturba- 
tively. To lowest order in a, the Green's functions read: 



Gu{e)^&']{ey 



- E 

n.n' —l,r,( 



(7) 



The self-energies are represented schematically in 

Fig. [1] by one-loop diagrams. The first diagram arises 



O 



E = 



+ 



+ 



FIG. 1: Diagrams for the self-energy. Dots denote vertices at 
tunneling from the single level to the leads, full and dashed 
lines stand for electron and phonon Green's functions, respec- 
tively. 

from the terms proportional to a^x'^ in the perturba- 
tive expansion of Eq. ^ with respect to a. It leads to 
a renormalization of tunneling matrix elements: ti^^ ^ 
ti/ril - a^(a;^)o/2) , where {x'^)o = coth[hLo/{2kBT)] for 
the unperturbed oscillator. The second diagram accounts 
for the shift of the oscillator's position due to the Lorentz 
force acting on it: {x}a = 4a/'-*'-'/(6^)- It does not con- 
tribute to the current to order o? . The third diagram 
contains the non-trivial part of the electron-phonon in- 
teraction. The sum of the three diagrams reads: 



.(2) 



E 



A ram - 



(8) 



where t,/^.d = i*,/^ = tii^{±ia{x) a- c? {x^)qI2), A^l, 

I, r) and zero otherwise. The components of an,n' (n, n' 
l,r,d) read 



GZ^'\t)D^{i) + G^:l^{t)D^'\t)_^ 
G^n^>)D''{t) + G^S^{t)D\t) 

■ (9) 



Here, D is the Green's function for unperturbed phonons 
(at a 0): D^/^{e) = 2hw/[{e ± ^0+)^ - {hwf] and 
D^{e) = -2iTT[8{e - ftiu) + S{e + huj)] coth[huj/{2kBT)]. 

Evaluation of the current up to terms by inserting 
eqs. (O-® into (|4]) is now straightforward. After lengthy 
calculations, we get / = /'^"' + a^/'^' -f . . . , where 



/(2) 



Gn 



(10) 



The coefficients 
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r°(£) =-2a(e) [a{e + huj) + a{e)] - 4j'^b{e)c{e,e + Hlu) - [26(e) - b{e + huj) - b{e - nuj)]nB 



(11a) 



(2(5,,,, - l)6(e)&(£ 
^{hwf -{e-e') 



e, e ■ 



hio) [b{e - huj) - b{e)] {hujf (.s,r, + s,,r,o /r^ 

2 {a{e) [b{e) + ^(e')] + 2s,,-fb[e)d[e, e')} , 



ntj) [6(e - fiij) + (1 - 25,,,0Ke)] 

(lib) 

(11c) 



are expressed through the functions a{e) = Re[rG^ 
b{e) = Im[rG^_,(e)], c{e,e') = a{e)a{e') + b{e)b{e'), 
and d(e,e') = a{e)b{e') — b{e)a{e'), while sj/^ ~ ili 
7 = (Fi — T]^)/T is an asymmetry factor for the cou- 
pling to the contacts, and ub = [e^^^''^'^ ~ is the 
Bose factor. The two first terms in the r.h.s. of Eq. (fTUl) 
express elastic as well as inelastic electron tunneling with 
emission/absorption of one phonon. The last term can- 
not be interpreted as a single-particle elementary process: 
It is related to the many-body character of the Fermi sea 
in the leads @, 0, @l- following, we discuss the 

result for the current in several regimes. 

Let us start with considering the linear conductance 
G = {dI/dV)y^^ = G(°) + a2G(2) + .... Far from the 
resonance, the tunnel regime is reahzed. Eqs. ©-([HI) 
yield (/i/e2)G^°l, = 4F,F,/£2 « l, where Sd ^Sd- {tM + 
lir)/'^- From Eq. (|10p . we obtain 



^tun ~ ^"-^tun 



1 



2nB - 2— —71 B[nB 
kRl 



1] 



(12) 



in agreement with Ref. That is, the magneto- 

conductance is negative at T = 0: AGtun = Gtun — 
G 



(0) 



-4a^Gt^'„, and it vanishes at high temperature 



like AGtun/G 



(0) 



-4a^ huj/{3kBT). 



The general expression of G*^^^ at T = is 



G(2 



-G 



ma x/ ^ 



sBsip) 



In 



2, ,2 



Y2 



1 



2s 



■ arctan 



+ ,(13) 



where A± (//) = a{n)a{iJL ± fiw) -|- 2j'^b{iJ,)c{ii, fi ± hw) and 
-S±(m) = o-il^)Kl^ i ^) + 2"f'^b{^j:)d{fi, ^ ± hhj). The cor- 
rection G^^^ is always negative. It results in a narrowing 
of the resonance line at finite magnetic field. The con- 
ductance peak at resonance (erf = 0) is: 



(7(2) _ _r< 



472r2 



■r2 + {hojy 



1 



2hjjj ( huj 

m — 

ttF V F 



(14) 



Remarkably the correction vanishes only for symmetric 
coupling between the dot and the leads (7 = 0), in con- 
trast with the polaronic problem where Gres = identi- 
cally. In the adiabatic limit, fku <^T, we find 

G(2) = -4G,naxr2[e2 ^ ^2p2]/[~2 ^ p2]2^ (^5) 
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FIG. 2: Gate-voltage dependence of the magneto- 
conductance at r = and different values of huj/T, for a 
nanotube contacted symmetrically (left) or asymmetrically 
(right). 



In the anti- adiabatic limit, fiw ^ T , the dominant term 
contributing to Eq. (fTS]) is 



G(2 



G,] 



-^212 



--d^- J • (16) 

This result receives a simple interpretation d: At high 
resonance frequency, the harmonic oscillator remains in 
its ground state, |0), and an effective Hamiltonian for 
electron tunneling in the device is obtained by project- 
ing Eq. ([T]) on this state. This describes a nonintcracting 
resonant level coupled to the leads through renormal- 
ized tunneling matrix elements t,(0|e*"''^'|0} = t^e~°' 
with corresponding level widths F,e~" . As a result, 
the resonance narrows, but the maximum transmission 
is unchanged. Eq. (|16p expresses this in lowest or- 
der in the coupling constant. The conductance reduc- 
tion arises to higher order in T/huj, AGres/Gmax oc 
~a^j'^T/{huj)ln{huj/T). The features discussed above 
arc clearly visible in Fig. [2l As in the polaronic trans- 
port problem Q , there is no vibrational sideband in the 
gate- voltage dependence of the linear conductance. 

It is worth mentioning that Eq. (fT3|) coincides (up 
to order a^) with a Landauer-Biittiker formula for the 
conductance in the elastic regime (zero bias-voltage 
and temperature) at arbitrary ratio hiu/T: G = 
{e^ /h)Y,kk'\Tki.k'r?5{£k)5{£k')- Here, the element of 
the T-matrix Tku^'r = (fffi)"'Gf; ^^^(ff^J-i is related 
to Gf; j.,^ that can be evaluated along the same lines as 
the Green's functions introduced above. We believe that 
this result is not so trivial. Indeed, contrarily to the 
polaronic Hamiltonian, the one that we are considering 
does not satisfy the "proportionate coupling" conditions 
to the right and left leads under which such a relation be- 
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tween conductance and transmission was proven in the 
presence of interactions Q . 

Does the effect discussed above really depend on the 
quantum vibrations of the oscillator? ^From Eq. PT?)) . 
we find that the magnetoconductance is suppressed when 
ksT 3> huj: 



AG{T) 

G(o)(r) 



3 si+r^ 



r < keT. 



This result differs from the one obtained with a classical 
description of the oscillator. In this case, x{t) obeys a 
classical equation of motion in presence of a Langevin 
force describing thermal fluctuations. The displacement 
x{t) acts as an external bias voltage T4c = 2a{h/e)x{t) 
in E g. dp . The current can then be calculated following 
Rcf. [10[. This classical contribution coincides with the 
term proportional to in the quantum result pOp - pi[l . 
It can be viewed as a rectification of Vac -fluctuations. We 
find that the corresponding magnetoconductance AGci 
is either negligible (for huj <C ksT <C F) or different 
(for huj,T <C ksT when it is of the same order) from 
the quantum result. In particular, in the tunnel regime, 
I{t) = Gtun[V + 2a{h/e)x{t)] and AGci = since (.t)ci 
vanishes. 

The formula pO|) for the current also allows to address 
the finite bias nonlinear regime. The possibility to excite 
the phonon (inelastic cotunncling) when eV > Huj leads 
to a nonanalytical voltage dependence of the current in 
vicinity of the threshold at T = 0, with leading terms: 



1 0/(2) _ ci F 



+ 4.020 {eV -hiu), (17) 



and asymptotic expressions for the coefficients 



^ -8fi^F3 (^^r. + si) I (F2 + sir , fi^ « F, 
^ -7(4F/M^£Vr, /!i^»F, 



and 



C2 



F2(72r2 + £2)/(r2 + ,~2)^^ ^<<r, 

(2F/M^(e2+72r2)/r2, fi^>F. 



A similar feature has been discussed for the polaron prob- 
lem It is clearly seen in the gate and bias voltage 
dependence of the magneto-conductance shown in Fig. [31 
In addition. Fig. [3] shows how the standard sequential 
tunneling lines at = ±ey/2 narrow at finite magnetic 
field, while additional phonon-assisted tunneling lines ap- 
pear at Ed = ±(ey/2 — hw). 

The weak-coupling approximation used above is well 
justified experimentally. However, it overlooks several 
features that may deserve further study. First, above 
the inelastic threshold [eV > hui), the phonon mode 
is driven out of equilibrium. By analogy with the po- 




FIG. 3: Gate and bias voltage dependence of the differential 
conductance dl^'^^ /dV at T = for Tiuj = lOF and F^. = F,. 



bias 5V <^ hw/e for fiij < F, and 5V < {hio) /Fe for 
hiij ^ F, under the additional condition that the current- 
induced damping of the oscillator min[(ftw)3/F2, F] re- 
mains much smaller than lo. Second, additional inelastic 
lines, corresponding to mu lti-p honon absorption, appear 
at id ~ ±{eV/2 — nhw) [I4I and lead to a magnetic- 
field-induced current suppression similar to the electro- 
static Frank-Condon blockade [l3| recently observed ex- 
perimentally [1^ . 

In conclusion, we have studied how the Aharonov- 
Bohm phase accumulated by the electrons crossing 
a vibrating nanotube affects its resonant magneto- 
conductance in the weak coupling regime. 
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